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Equalization  in  Aeronautical  Telemetry 
Using  Multiple  Antennas 

Michael  Rice,  Md.  Shah  Afran,  Mohammad  Saquib 


Abstract 

This  paper  introduces  a  generalized  time-reversed  space-time  block  codes  (GTR-TRSTBC)  as  both 
a  conceptual  tool  to  examine  the  impact  of  unequal  power  allocation  in  equalized  single-carrier  systems 
employing  two  transmit  antennas  and  as  a  method  for  exploiting  partial  channel  state  information  by  the 
transmitter.  The  generalization  involves  unequal  power  allocation  parameterized  by  p.  The  criteria  for 
selecting  the  optimum  p  is  minimizing  the  residual  mean-squared  error  at  the  MMSE  equalizer  output. 
GTR-STBC  is  applied  to  measured  channel  impulse  responses  and  a  simple  statistical  channel  model. 
The  results  show  that  1)  the  optimum  value  of  p  gives  the  best  tradeoff  between  signal-to-noise  ratio  and 
ISI;  2)  equal  power  allocation  may  not  be  the  optimum  power  allocation  when  channel  state  information 
is  available;  and  3)  the  optimum  profile  of  p  over  measured  channels  is  significantly  different  than  from 
that  in  statistical  channel  models. 


I.  Introduction 

Aeronautical  telemetry  is  a  wideband  air-to-ground  communication  link  with  severe  size, 
weight  and  power  limitations  imposed  on  the  airborne  transmitter  combined  with  frequency 
selective  fading  channels.  The  size,  weight,  and  power  constraints  on  the  airborne  platform 
impose  limitations  on  the  transmitter.  Many  of  these  needs  are  met  by  using  highly  power- 
efficient  transmitters,  such  as  radio  frequency  (RF)  power  amplifiers  operating  in  full  saturation. 
Consequently,  given  constant-envelope  modulations  are  preferred.  This  preference  manifests  itself 
by  the  modulations  described  in  the  aeronautical  telemetry  standard  IRIG  106  [1]:  pulse  code 
modulation  with  frequency  modulation  (PCM/FM),  shaped  offset  quaternary  phase  shift  keying 
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version  TG  (SOQPSK-TG)  and  a  version  of  non-binary,  multi-index  continuous  phase  modulation 
(CPM)  called  ARTM  CPM.  PCM/FM  is  the  legacy  modulation  which,  relative  to  the  other  two,  is 
not  as  bandwidth  efficient.  The  importance  of  SOQPSK-TG  has  increased  over  the  past  15  years 
as  the  availability  of  bandwidth  has  come  to  dominate  scheduling  and  regulatory  discussions. 

Multipath  fading  tends  to  be  the  dominate  impairment  in  maintaining  link  availability.  Given 
the  high  data  rates  required  by  modem  test  articles,  the  corresponding  RF  bandwidth  of  the 
telemetry  signal  is  such  that  essentially  all  multipath  fading  is  frequency  selective.  Because 
frequency  selective  fading  leads  to  inter-symbol  interference  (ISI),  techniques  for  mitigating  ISI 
using  the  IRIG  106  modulations  (SOQPSK-TG  in  particular)  are  of  immediate  interest. 

The  goal  therefor  is  one  involving  single-carrier,  constant-envelope  modulation  operating  in  a 
frequency  selective  channel.  This  is  clearly  an  equalization  and/or  diversity  problem.  Diversity 
reception,  especially  by  widely  separated  antennas,  is  an  obvious  solution.  But  given  the  expense 
of  ground-based  receiving  installations,  this  solution  has  been  relatively  unpopular  at  the  major 
test  range  facilities.  On  the  other  hand,  transmit  diversity  can  be  a  workable  solution  when 
size,  weight,  and  power  constraints  allow  it.  The  most  common  scenario  will  be  one  involving 
multiple  transmit  antennas  with  a  single  receive  antenna. 

It  is  well  known  that  multiple  antenna  systems  are  capable  of  increasing  reliability  or  through¬ 
put  in  multipath  fading  channels.  In  flat  fading,  the  optimum  signaling  approach  depends  on 
what  the  transmitter  knows  [2].  If  the  transmitter  knows  the  channels  between  each  transmit 
and  receive  antenna,  then  spatio-temporal  coding  [3]  is  optimum  in  that  it  maximizes  signal- 
to-noise  ratio  [2].  If  the  transmitter  does  not  have  this  knowledge,  then  a  diversity-maximizing 
orthogonal  design  (such  as  the  Alamouti  code  [4])  is  optimum  [2].  In  frequency  selective  fading, 
the  general  approach  is  to  use  OFDM  and  apply  these  techniques  on  a  per  subcarrier  basis. 
Given  the  constraints  described  above,  OFDM  is  not  a  viable  approach  and  approaches  suitable 
to  single-carrier  modulation  are  of  interest. 

In  the  case  of  multi-transmit-antenna  systems  employing  single-carrier  modulation  and  operat¬ 
ing  on  frequency  selective  channels,  time-reversed  space-time  block  codes  (TR-STBCs)  [5]-[9] 
play  the  role  that  the  Alamouti  code  does  in  frequency  non-selective  fading  with  an  uninformed 
transmitter.  When  the  transmitter  knows  the  channels,  the  situation  becomes  less  clear.  Given  the 
fact  that  SOQPSK-TG  does  not  have  any  “variable  parameters”  such  as  number  and  location  of 
constellation  points,  just  how  an  informed  transmitter  is  to  use  channel  state  information  is  not 
obvious.  The  only  other  variable  available  to  the  system  is  the  power  allocated  to  each  transmit 
antenna. 

In  an  effort  to  explore  the  possibility  of  transmit  diversity  to  increase  the  robustness  of 
aeronautical  telemetry  links  prone  to  multipath  fading,  channel  sounding  experiments  involving 
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multiple  transmit  antennas  have  been  eondueted  at  Cairns  Army  Airfield,  Ft.  Rueker,  Alabama, 
[10],  [11]  and  Edwards  Air  Foree  Base,  California,  [20].  Initial  experiments  involving  TR-STBC 
applied  to  the  ehannel  the  helieopter-to-ground  ehannel  sets  revealed  some  eurious  behavior  This 
behavior  with  TR-STBC,  QPSK,  and  the  ehannel  sets  was  deseribed  in  [11].  This  same  behavior 
TR-STBC,  SOQPSK-TG  for  the  same  data  set  is  illustrated  in  Fig.  1.  This  plot  eompares  the 
simulated  bit  error  rate  (BER)  performanee  of  TR-STBC  using  MMSE  equalizers  (deseribed  in 
Seetion  III)  to  the  simulated  BER  performanee  of  a  link  using  only  one  of  the  two  available 
ehannels.  TR-STBC  performs  better  than  the  single-ehannel  link  using  only  ehannel  2,  but  worse 
than  the  single-ehannel  link  using  only  ehannel  1.  Clearly,  ehannel  1  is  better,  in  some  sense, 
than  ehannel  2.  In  faet,  ehannel  1  is  so  mueh  better  than  ehannel  2  that  ineorporating  ehannel 
2  into  a  TR-STBC  system  only  makes  things  worse. 

In  other  words,  there  are  cases  where  it  is  better  to  abandon  traditional  TR-STBC  and  use  only 
one  of  the  two  available  channels.  The  fact  that  this  curious  behavior  can  occur  on  real  channels 
prompts  one  to  ask,  “Can  this  curious  behavior  be  derived  from  the  given  impulse  responses  of 
two  channels?”  In  other  words,  it  appears  that  with  some  form  of  channel  state  information,  it 
is  possible  to  achieve  better  performance  than  TR-STBC. 

This  paper  answers  the  question.  As  a  conceptual  tool,  we  consider  a  fixed-power  trans¬ 
mitter  that  allocates  a  portion  of  the  fixed  power  to  each  channel.  This  power  allocation  is 
parameterized  by  p,  the  proportion  of  total  power  allocated  to  channel  1.  Using  only  one  of 
the  two  available  channels  is  captured  by  the  case  p  =  1  (channel  1  only)  or  p  =  0  (channel 
2  only).  Next,  we  introduce  Generalized  TR-STBC  (GTR-STBC),  a  modified  version  of  TR- 
STBC  that  incorporates  the  unequal  power  allocation  p.  Because  TR-STBC  and  its  generalized 
version  involve  the  combination  of  diversity  processing  and  equalization,  a  suitable  equalization 
technique  for  SOQPSK-TG  is  needed.  We  develop  and  analyze  an  minimum  mean-squared  error 
(MMSE)  equalizer  that  operates  on  the  samples  of  SOQPSK-TG.  Equipped  with  GTR-STBC 
and  a  suitable  equalization  technique,  we  use  the  mean-squared  error  criterion  to  identify  the 
optimum  p.  The  motivations  for  the  mean-squared  error  criterion  are  the  following: 

1)  Mathematically  tractable  expressions  for  the  mean-squared  error  at  the  equalizer  output 
are  easily  derived  [see  (42)  below]. 

2)  The  mean-square  error  criterion  neatly  captures  the  contributions  of  both  ISI  and  additive 
noise  at  the  equalizer  output. 

3)  Generally  speaking,  reductions  in  mean-squared  error  lead  to  reductions  in  bit  error  rate. 
By  expressing  the  residual  mean-squared  error  at  the  equalizer  output  as  a  function  of  p,  we  are 
able  to  chose  p  to  minimize  the  residual  mean-squared  error  for  a  given  pair  of  channel  impulse 
responses. 
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Fig.  1.  Simulated  BER  plots  for  10-Mbit/s  SOQPSK-TG  using  MMSE  equalizers  over  a  pair  of  representative  impulse 
responses  measured  over  a  helicopter-to-ground  channel.  The  circular  markers  are  the  BER  results  using  a  single  antenna,  the 
square  markers  are  the  BER  results  using  hoth  channels  with  TR-STBC. 


Finally,  we  apply  the  eoncepts  to  the  set  of  measured  channel  impulse  responses  to  see  if 
the  mean-squared  error  criterion  is  capable  of  identifying  the  cases  where  the  curious  behavior 
occurs.  The  results  show  that  the  mean-squared  error  criterion  does  indeed  capture  the  curious 
behavior. 

The  conceptual  tool  of  unequal  power  allocation  together  with  GTR-STBC  define  a  simple 
transmit  diversity  scheme  based  on  the  partial  knowledge  of  the  channel  by  the  transmitter.  Here, 
the  transmitter  only  needs  to  know  p,  which  is  easy  to  compute  at  the  receiver  and  send  back 
to  the  transmitter.  This  simple  scheme  includes  transmit  selection  diversity  (p  =  0  or  1)  and 
traditional  TR-STBC  (p  =  ^/2). 

This  paper  is  organized  as  follows.  SOQPSK-TG  is  described  in  Section  II  and  the  MMSE 
equalizer  for  SOQPSK-TG  is  developed  in  Section  III.  GTR-STBC  is  introduced  and  analyzed 
in  Section  IV.  GTR-STBC  is  applied  to  the  measured  channel  impulse  responses  in  Section  V, 
and  conclusions  are  presented  in  Section  VI. 
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II.  SOQPSK-TG 

Shaped  offset  QPSK,  version  TG  (SOQPSK-TG)  is  constant-envelope  modulation  with  good 
spectral  containment  and  good  detection  efficiency.  For  these  reasons,  it  is  an  ideal  choice  for 
aeronautical  telemetry  and  has  been  a  part  of  the  IRIG-106  standard  since  2004  [1].  SOQPSK-TG 
is  a  continuous  phase  modulation  (CPM)  whose  complex-valued  baseband  equivalent  represen¬ 
tation  [12]  is 

f2^E~ 

x{t,(x)  =  ^^exp[j{(j){t,(x)  +  (j)o)].  (1) 

The  phase  is 


0(f,  ck)  =  2nh 


OO 

an  fir  -  nTb)dT 

n=—oo 


OO 

=  2nh  anq(t  —  nTb). 

n=—oo 


(2) 


where  f{t)  is  the  frequency  pulse;  q{t)  =  f{'T)dT  is  the  phase  pulse;  0o  is  an  arbitrary  phase 
which,  without  loss  of  generality,  can  be  set  to  0;  is  the  bit  interval  (or  reciprocal  of  the  bit 
rate);  h  =  ^/2  is  the  modulation  index;  and  an  €  {  —  1,0, 1}  are  the  ternary  symbols  which  are 
related  to  the  binary  input  symbols  G  (  —  1, 1}  by  [13] 

=  {an-  an-2)  ^  ^3^ 

The  frequency  pulse  for  SOQPSK-TG  is  described  in^  [15],  [16].  Because  the  frequency  pulse 
spans  8  bit  times,  SOQPSK-TG  is  a  partial  response  CPM.  Consequently,  the  maximum  likeli¬ 
hood  detector  in  the  additive  white  Gaussian  noise  (AWGN)  environment  is  a  sequence  estimator 
[17],  usually  implemented  using  the  Viterbi  algorithm  involving  a  trellis  with  512  states.  A 
reduced-complexity,  near  optimum  4-state  detector  was  described  in  [15].  The  most  common 
detector  applied  to  SOQPSK-TG  is  the  simple  OQPSK-type  symbol-by-symbol  (SxS)  detector 
described  by  Perrins  [18].  The  SxS  detector  comprises  a  detection  filter  (see  Fig.  3  of  [18]) 
whose  output  is  sampled  at  1  sample/bit  and  is  treated  as  offset  QPSK. 


III.  MMSE  Equalization  of  SOQPSK-TG 

The  system  considered  here  is  summarized  in  Eigure  2.  As  before,  the  complex-valued  base¬ 
band  equivalent  representation  [12]  is  used  for  all  signals.  Starting  with  the  block  diagram  of 


*The  basic  form  was  first  introduced  by  Hill  in  [14]  where  “-A”  and  “-B”  versions  were  described.  SOQPSK-A  has  a  slightly 
narrower  bandwidth  (measured  at  the  -60  dB  level)  and  slightly  worse  detection  efficiency  than  SOQPSK-B.  The  Telemetry 
Group  (TG)  of  the  Range  Commanders  Council  adopted  the  compromise  waveform,  designated  SOQPSK-TG  in  2004. 
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Figure  2  (a),  the  SOQPSK-TG  signal  Xc(t)  is  transmitted  through  a  channel  impulse  response 
hc{t)  whose  output,  accompanied  by  thermal  noise,  forms  the  receive  signal  rdt).  After  the 
application  of  an  anti-aliasing  low-pass  filter  with  impulse  response  ha{t),  T-spaced  samples  of 
rc{t)  are  produced  by  an  A/D  converter.  Assuming  the  anti-aliasing  filter  does  not  distort  the 
received  signal,  the  samples  of  the  received  signal  may  be  expressed  as 

L2 

r{n)  =  x{n)  *  h{n)  -f  w{n)  =  h{k)x{n  —  k)  +  w{n)  (4) 

k=-Li 


where 

r{n)  =  r-ciriT)  x{n)  =  Xc{nT)  h{n)  =  hc{t)  * 


and  where  the  w{n)  is  the  n-th  sample  in  a  sequence  zero-mean  complex- valued  Gaussian 
random  variables  with  autocorrelation  function 

Rw{k)  =  ^^w{n)w*{n  -  k)^  =  (Tl6{k).  (5) 

Note  that  (4)  assumes  the  discrete-time  channel  has  support  on  —Li  <n<  L2.  As  a  first  step, 
we  apply  the  samples  r(n)  to  a  filter  matched  to  the  disccrete-time  channel  to  produce  y{n): 


(6) 

(7) 


where  Lgq  =  Li  +  L2  and  v{n)  is  a  complex  valued  Gaussian  random  sequence  with  zero  mean 
and  autocorrelation  function 


The  samples  y{n)  form  the  input  to  an  MMSE  equalizer.  The  MMSE  equalizer  is  an  EIR  filter 
with  coefficients  c(n)  for  —Lc  <n<Lc  designed  to  minimize  the  mean  squared  error  between 
the  equalizer  filter  output  x{n)  and  the  sequence  x{n).  The  entire  system  may  be  represented 
by  an  equivalent  discrete-time  system  shown  in  Eigure  2  (b). 

The  challenge  with  equalizing  samples  of  the  modulated  signal  is  that  the  underlying  continuous¬ 
time  waveform  is  not  wide-sense  stationary  [12].  This  fact  carries  over  to  the  samples  of  Xc{t) 
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(9) 


and  is  captured  by  the  fact  that  the  autocorrelation  function  of  the  x{n)  is  of  the  form 

that  is,  the  autocorrelation  function  is  a  function  of  both  sample  indexes,  not  the  difference 
between  them.  Consequently,  the  equalizer  filter  coefficients  are  a  function  of  the  sample  index 
n.  It  is  hard  to  see  how  this  solution  has  any  practical  utility,  especially  in  the  presence  of  a 
real-time  performance  requirement.  In  the  end,  the  designer  is  left  with  suboptimal  approaches 
of  reduced  computational  complexity  whose  accompanying  performance  penalty  is  acceptable. 

The  simplest  suboptimal  approach  is  to  assume  the  signal  samples  are  wide-sense  stationary. 
Here,  the  autocorrelation  function  is  of  the  form 

RAk-i)  =  ^E{x{k)x*{i)},  (10) 

that  is,  the  autocorrelation  function  depends  on  the  difference  of  the  sample  time  indexes.  The 
wide-sense  stationary  assumption  for  x{n)  greatly  simplifies  the  solution.  Because  the  optimum 
equalizer  coefficients  no  longer  depend  on  the  samples  index  n,  the  relationship  between  x{n) 
and  the  equalizer  output  x{n)  is 

x{n)  =  c{n)  *  y{n)  =  c{k)y{n  —  k).  (11) 

/c=  Lq 


The  vector  of  filter  coefficients  that  minimizes  the  mean  squared  error 


S  =  E 


—  x{n) 


2 


(12) 


is  given  by 


c  = 


-f  R^, 


Rx,2g'^ 


(13) 


where  c  is  the  (2Lc  -f- 1)  x  1  vector  of  filter  coefficients,  G  is  the  (2Lc  -f  1)  x  (2Lc  -f  1  -f  2Leq) 
matrix 

r^eq(f^c)  ■  ■  ■  ^eq(  1 


G  = 


^eq(f^c)  ■  ■  ■  ^eq( 


(14) 


^eq(f^c)  ■  ■  ■  ^eq( 
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R^.i  is  the  (2Lc  +  2Leq  +  1)  x  (2Lc  +  2Leq  +  1)  matrix 


Rx,i 


Rx(0) 

Rxil) 


Rx(-l) 

Rx(0) 


RJ-2Lr  -  2L, 


eqy 


Rx{~‘^Lc  —  2Leq  +  1) 


-Rx(2i^c  +  2Leq)  Rx{2Lc  +  2Lf.^  —  1)  ■■■  -Rx(O) 

Rx  is  the  (2Lc  +  1)  x  (2Lc  +  1)  noise  autocorrelation  matrix  given  by 

^  R„(0)  ■  ■  ■  R„(-2L,  ^ 


Rt  = 

[Rv{2L, 

Rx,2  is  the  (2Lc  +  1)  x  (2Lc  +  1)  matrix  given  by 

■  R,(0)  R,(-l) 

-Rx(l)  -Rx(O) 

Rx,2  = 

Rx{‘^Lc)  Rx{‘^Lc  +  1) 

and  g  is  the  1  x  (2Lc  +  2Leq  +  1)  vector  given  by 


Rv{0)  J 


-Rx(~2Lc) 

Rx{~‘^Lc  +  1) 

-Rx(O) 


g  = 


Olx(Lc-Leq)  heq  Olx(L,-L,q) 


(15) 


(16) 


(17) 


(18) 


where  Oix{Lc-l^)  is  a  row  vector  comprising  Lc  —  L^q  zeros  (here  we  assume  Lc  >  L^q,  i.e.,  the 
equalizer  is  longer  than  the  channel),  and  hgq  is  the  1  x  (2Leq  +  1)  vector  given  by 


^eq 


hsqi^L, 


eq; 


/ieq(O) 


^eq(  -^eq) 


(19) 


The  question  is  now,  what  function  should  be  used  for  the  autocorrelation  function  Rx{k)2  Two 
approximations  are  investigated  here.  The  first  is  an  empirically-derived  autocorrelation  function. 
The  empirical  autocorrelation  function  is  obtained  by  generating  a  large  number  of  samples  x{n) 
and  using  the  standard  estimation  technique  assuming  wide  sense  stationarity.  Given  L  samples 
of  a;(n)  for  n  =  0, 1, . . . ,  L  —  1,  this  empirical  autocorrelation  function  is 


R,{k) 


1 

2{L-k) 


L-l 

x{n)x*{n  —  k), 

n=k 


0<k<L-l 


(20) 


together  with 


R,{k)=Rl{-k), 


-L  <  k  <0. 


(21) 
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delay 


Fig.  3.  A  plot  of  the  empirical  autocorrelation  function  for  SOQPSK-TG:  (top)  the  real  part  of  i?c(fc);  (hottom)  the  imaginary 
part  of  7?c(fc).  The  sample  rate  for  the  SOQPSK-TG  samples  is  at  2  samples/hit.  Markers  indicate  the  values  for  —5  <  k  <  5. 


A  plot  of  Re{k)  corresponding  to  L  =  2  x  10®  samples  of  SOQPSK-TG  sampled  at  2  samples/bit 
is  shown  in  Figure  3  for  the  first  100  lags  (i.e.,  —100  <  k  <  100).  The  top  plot  shows  the  real 
part  of  Rs{k)  and  the  lower  plot  shows  the  imaginary  part  of  Re{k).  The  only  signifieant  values 
are  those  for  —5  <  A;  <  5  and  indieated  by  markers  on  the  plot.  Consequently,  in  the  simulation 
results  presented  below,  we  assume  Re{k)  =  0  for  \k\  >5. 

The  seeond  approximation  is  to  assume  the  data  are  uneorrelated.  This  generates  a  eorrelation 
funetion  of  the  form 


Ri{k)  =  al5{k). 


(22) 
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Here,  the  eorresponding  eorrelation  matriees  and  R^,  2  are  sealed  identity  matriees  and  they 
funetion  as  regularizers  in  the  numerieal  eomputations  (13). 

To  test  the  performanee  of  these  approximations,  we  simulate  the  performanee  of  the  equalized 
SOQPSK-TG  system  outlined  in  Figure  4.  For  this  test,  we  use  a  ehannel  obtained  from  our 
ehannel  measurement  eampaign  deseribed  below  with  Lgq  =  19  and  operating  at  iV  =  2 
samples/bit.  The  ehannel  transfer  funetion  is  plotted  in  Figure  5.  The  equalizer  output  is  applied 
to  the  symbol-by-symbol  deteetor  for  SOQPSK-TG  (denoted  “SxS  deteetor”  in  the  figure) 
deseribed  in  Seetion  IT 

The  effeetiveness  of  using  Re{k)  and  Ri{k)  in  eomputing  the  equalizer  eoeffieients  was 
measured  by  the  simulated  bit  error  rate.  The  simulations  applied  the  equalizer  eoeffieients 
obtained  from  (13)  for  Lg  =  5  x  L^q  =  95  [19]  and  using  Rs{k)  and  Ri{k).  The  simulated 
bit  error  rate  is  plotted  in  Figure  6,  where  we  observe  that  there  is  essentially  no  differenee  in 
the  performanee  between  the  two.  Consequently,  we  prefer  Ri{k)  over  Re{k)  in  what  follows 
beeause  this  ehoiee  simplifies  the  eomputations  of  the  equalizer  filter  eoeffieients.  In  this  case, 
the  optimum  equalizer  filter  coefficients  are  given  by 

where  Hgq  is  the  {2Lc  -f-  1)  x  (2Lc  -(- 1)  matrix  given  by 

■  /igq(O)  ■■■  /lgq(-2Lc)' 

Hgq  =  ;  ;  (24) 

_/lgq(2Lc)  ■  ■  ■  /igq(O) 

where  it  is  understood  that  hsq{k)  =  0  for  |/c|  >  Lgq.  The  corresponding  mean  squared  error  is 

^min  =  Cr^  1^1  -  g  GGt  +  ^^^^Hgq^  gt ^  .  (25) 
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2  samples/bit 


1  sample/bit 


Fig.  4.  A  block  diagram  of  the  system  used  to  test  the  performance  of  the  SOQPSK-TG  autocorrelation  assumptions. 


Fig.  5.  The  channel  used  in  the  simulations  reported  in  Figure  6.  The  channel  is  obtained  from  channel  measurements  described 
below.  The  bit  rate  is  10  Mbits/s  and  the  channel  is  sampled  at  2  samples/bit.  The  corresponding  feeq(n)  has  Leq  =  19- 
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E^/N^(dB) 


Fig.  6.  Simulated  bit  error  rate  results  for  the  system  shown  in  Figure  4  using  the  channel  illustrated  in  Figure  5.  Also  included 
for  reference  is  the  performance  of  the  symhol-hy-symhol  detector  in  the  AWGN  channel  and  the  optimum  maximum  likelihood 
sequence  detector  (MLSD)  in  the  AWGN  channel  (see  [18]). 
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IV.  Generalized  TR-STBC  (GTR-STBC):  Non-Equal  Power  Allocations  Using 

TR-STBC 

An  abstraction  (to  the  sample  level)  for  a  2  x  1  TR-STBC  system  is  illustrated  in  Fig.  7. 
Here  the  system  transmits  2L  samples  of  an  SOQPSK-TG  signal  sampled  at  N  samples/bit  over 
two  transmit  antennas  to  one  receive  antenna.  The  equivalent  discrete-time  channel  between 
transmit  antenna  1  and  the  receive  antenna  is  represented  by  the  impulse  response  hi{n)  for 
—-^1,1  <  ^  <  Li  2  whereas  the  equivalent  discrete-time  channel  between  transmit  antenna  2  and 
the  receive  antenna  is  represented  by  the  impulse  response  h2{n)  for  — L2,i  <  n  <  L2,2- 

The  TR-STBC  encoder  partitions  the  sample  sequence  a;(0), . . . ,  x{2L  —  l)  into  two  sequences 
Xi{n)  and  X2{n)  as  shown  in  Fig.  7.  The  length-2L  packet  is  transmitted  in  two  intervals^  each 
spanning  L  sample  intervals.  During  the  first  interval  a;i(0), . . . ,  a;i(L  —  1)  is  transmitted  from 
antenna  1  whereas  0:2(0), ...,  0:2 (T  —  1)  is  transmitted  from  antenna  2.  During  the  second  interval, 
xl{L  —  l), ...,  0:2(0)  is  transmitted  from  antenna  1  whereas  —xl{L—l), . . . ,  —0:^(0)  is  transmitted 
from  antenna  2. 

Power  division  using  0  <  p  <  1  is  accomplished  by  modified  the  TR-STBC  system  along  the 
lines  illustrated  in  Fig.  7.  Amplitude  scaling  is  applied  to  the  signals  entering  each  channel  so 
as  to  divide  the  power  between  the  channels.  Here  p  represents  the  proportion  of  total  power 
allocated  to  transmit  antenna  1.  The  traditional  TR-STBC  system^  is  a  special  case  for  which 
p  =  y2.  The  square-root  is  used  in  Fig.  7  because  the  amplitudes  are  what  are  being  modified — 
the  energy  (or  power)  is  the  square  of  the  amplitude. 

The  received  signal  r(n)  is  given  by 

r{n)  =  ^/pXl{n)  *  hi{n)  +  a/1  -  px2{n)  *  h2{n)  +  w{n)  (26) 


where  w{n)  is  a  complex- valued  Gaussian  random  sequence  with  zero  mean  and  autocovariance 
function  given  by  (5).  The  TR-STBC  decoder  partitions  r(n)  into  ri(n)  and  r2(n)  as  follows: 


ri(n)  =  r(n)  for  0  <  n  <  L  —  1 
r2(n  —  L)  =  r{n)  for  L  <  n  <  2L  —  1. 


(27) 


^In  a  practical  implementation,  a  guard  interval  at  least  as  long  as  the  longest  channel  impulse  response  must  be  inserted 
between  the  two  intervals.  Here,  such  an  interval  is  assumed,  although  we  won’t  complicate  the  notation  to  make  this  explicit. 

^In  the  traditional  TR-STBC  system,  p  =  1/2  is  included  in  neither  the  development  nor  the  notation  because  the  same  power 
is  assumed  to  be  applied  to  each  channel.  Hence  there  is  no  need  to  account  for  it,  other  than  in  normalizing  the  noise  variance. 
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These  two  sequences  are  given  by 

^i(^)  =  a/p  xi{n)  *  hi{n)  +  a/1  -  p  X2{n)  *  /i2(n)  +  wi{n)  (28) 

r2{n)  =  ^p  xl{-n)  *  hi{n)  -  a/1  -  p  *  h2{n)  +  (29) 

where  wi{n)  and  W2{n)  are  related  to  w{n)  the  same  way  ri(n)  and  r2(n)  are  related  to  r(n). 
The  TR-STBC  decoder  processes  ri(n)  and  r2(n)  using  a  bank  of  filters  based  on  the  channel 

impulse  responses  hi{n)  and  h2{n)  as  shown.  The  result  of  this  processing  is  a  pair  of  parallel 

sequences  yi{n)  and  y2{n)  which  may  be  expressed  as 


pi(n)  =  ri(n)  *  y/p  hl{-n)  -  r^i-n)  *  /i2(n) 

=  Xi{n)  *  [p  hi{n)  *  hl{—n)  +  (1  —  p)  h2{n)  *  h*Y—n) 


hcq{n) 


+  Wi{n)  *  y/p  hl{-n)  +  W2{-n)  *  a/1  -  p  h2{n)  (30) 


vi{n) 


y2{n)  =  ri(n)  *  ^/l-  p  hl{-n)  +  r*2{-n)  *  y/^  hi{n) 

=  X2{n)  *  [(1  -  p)  /i2(n)  *  hl{-n)  +  p  hl{-n)  *  hi{n) 


hcq{n) 


+  a/1  -  p  wi{n)  *  h2{-n)  +  y/p  W2{-n)  *  hi{n) .  (31) 

" - V - ' 

V2{n) 

These  equations  show  that  the  equivalent  composite  channel  for  non-equal  power  allocation  is 


heq{n)  =  p  hi{n)  *  hl{-n)  +(1  -  p)  /i2(n)  *  hl{-n) .  (32) 

^ - V - '  ^ - V - ' 

rii{n)  V2{n) 

Because  the  support  for  hi{n)  is  — Ti  i  <  n  <  Li  2,  the  support  for  pi(n)  is  —  (Li  i+Li  2)  <  n  < 
(Li  1  +  Li  2).  Similarly  because  the  support  for  h2{n)  is  —^2,1  <n<  ^2,2,  the  support  for  p2(^) 
is  —(^2,1  +  ^2,2)  <n<  (L2,i  +  ^2,2)-  Consequently,  the  support  for  /leq(rt)  is  — Lgq  <n<  L^q 
where 

Ceq  =  max  I  Ti  l  +  Ti^2)  T2,l  +  ^2,2  I  •  (33) 

The  noise  sequences  vi{n)  and  V2{n)  are  complex- valued  Gaussian  random  sequences  each  with 
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zero  mean  and  autocorrelation  and  cross  correlation  functions 

l^vi{n)vl{n  -  A;)|  =  |t;2(n)t;;(n  -  /i:)|  =  alK^{k).  (34) 

E  |t;i(n)w2(ra  — /c)|  =  0.  (35) 

By  way  of  summary,  the  TR-STBC  system  presents  to  the  equalizers  the  sequences  yi  (n)  and 
7/2(77)  which  may  be  represented  by 

7/1(77)  =  0:1(77)  *  heq{n)  +  7:1(77)  (36) 

7/2 (n)  =  0:2(77)  *  heq{n)  +  7:2(77)  (37) 

where  heq{n)  is  given  by  (32).  The  noise  terms  7:1(77)  and  7:2(77)  are  uncorrelated  zero-mean 
Gaussian  random  sequences  each  with  autocorrelation  function  (34).  A  pair  of  equalizers  operate 
in  parallel  on  7/1(77)  and  7/2  (^^)-  Because  the  the  noise  sequences  7:1(77)  and  7:2(77)  are  statistically 
equivalent  and  hsq{n)  is  common  to  both,  the  pair  of  equalizers  operating  on  7/1(77)  and  7/2(77)  are 
identical  as  long  as  01(77)  and  02(77)  are  statistically  equivalent  (the  usual  case).  Any  equalizer 
can  be  applied  here  (linear  or  non-linear,  with  or  without  noise  whitening)  with  the  usual 
performance-complexity  tradeoffs.  Here,  we  apply  the  approximate  MMSE  equalizer  introduced 
in  Section  III  because  the  MMSE  equalizer  permit  a  mathematically  tractable  analysis  for  the 
resulting  mean-squared  error.  We  leverage  the  analytical  expression  to  find  the  value  of  p  that 
minimizes  the  mean  squared  error. 

The  equalizer  filters  of  Eig.  7  are  identical,  and  the  vector  of  filter  coefficients  is  given  by  (23) 
and  the  corresponding  mean  squared  error  is  given  by  (25).  In  these  equations,  G,  g,  and  Hgq 
may  be  expressed  in  terms  of  power  allocation  parameter  p.  Using  (32),  it  is  straightforward  to 
show  that 


G  =  pGi  +  (l-p)G2  (38) 

g  =  pgi  +  (1  -  p)g2  (39) 

Heq=pHi  +  (l-p)H2  (40) 

where  Gi,  gi,  and  Hi  are  formed  from  771(77)  the  same  way  G,  g,  and  Hgq  are  formed 
from  hsq{n),  respectively.  Similar  definitions  apply  to  G2,  g2,  and  H2  with  772(77).  Making 
the  substitutions  for  G,  g,  and  Hgq  gives 

c  =  M-i  (^pgi  +  (1  -p)g2y 


(41) 
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and 


C  _  — i 

^min  ^  X 


1-  pgl  +  (l-p)g2  M  pgl  +  (l-p)g2 


(42) 


where 


M=  (pGi  +  (l-rtG2)  (pGi  +  (l-rtG2)’  +  ^l^^(pHi  +  (l-p)H2),  (43) 


Equation  (42)  is  the  desired  relationship:  for  a  given  pair  of  channels  hi{n)  and  /i2(n)  ,  it  expresses 
the  mean  squared  error  at  the  output  of  the  MMSE  equalizer  as  a  function  of  the  power  allocation 
p.  Thus,  for  a  fixed  pair  channels,  one  can  choose  the  power  allocation  to  minimize  the  achievable 
mean-squared  error. 


V.  Numerical  Results 


The  forgoing  analysis  was  applied  to  measured  impulse  responses  from  four  multipath  en¬ 
vironments  typically  encountered  in  aeronautical  telemetry.  Eor  the  modulation,  we  assume  10 
Mbit/s  SOQPSK-TG  with  the  equalizers  operating  at  2  samples/bit.  In  each  of  the  four  cases, 
we  compute  the  optimum  power  allocation  parameter  p  to  a  pair  of  channels  labeled  hi  (n)  and 
h-2{n).  The  physical  meaning  of  these  channels  depends  on  the  scenario  as  described  below. 

Two  normalizations  are  applied  to  the  channels:  the  natural  normalization  and  the  equal- 
energy  normalization.  Let  /ii,u(rt)  and  h2^a{n)  be  unnormalized  channel  impulse  responses  for 
the  two  equivalent  discrete-time  channels  obtained  directly  from  the  channel  sounding  data,  and 
let 

Ti.2  T2.2 

Ei=  Y.  E2=  Y  l^2,u(n)P  (44) 

n=-Li_i  n=— L2,i 


be  the  energies  in  two  channels.  The  natural  normalization  uses 


hi{n)  =  -^/ii,u(n)  h2{n)  =  -^/i2,u(n) 


(45) 


where  E  =  max{Ei,  E2}.  This  normalizes  the  stronger  of  the  two  channels  to  unit  energy."^  We 
call  this  the  natural  normalization  because  in  real  multi-antenna  scenarios,  it  is  often  the  case 
that  one  of  the  channels  is  stronger  than  the  other. 

Eor  the  equal-energy  normalization,  we  use 


hi{n)  = 


h2{n) 


1 

\/E2 


h2,viin). 


''The  motivation  for  unit  energy  is  for  scaling  the  noise  variance  to  define  the  signal-to-noise  ratio. 


(46) 
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Here,  both  ehannels  are  normalized  to  unit  energy.  This  is  more  typieal  of  statistieal  or  math- 
ematieal  models  of  multi-antenna  propagation  and  is  a  eommon  feature  of  the  more  aeademie 
approaeh. 

The  numerieal  results  were  produeed  as  follows.  For  eaeh  pair  of  ehannel  impulse  responses, 
the  impulse  responses  were  normalized  using  one  of  the  two  proeedures  deseribed  above.  The 
value  of  p  that  minimizes  (42)  for  Lc  =  5  x  Leq  [19]  was  eomputed. 

A.  Cairns  Army  Airfield,  Ft.  Rucker,  Alabama 

The  first  data  set  eomprises  ehannel  impulse  responses  eolleeted  at  Cairns  Army  Airfield, 
Ft.  Rueker,  Alabama.  Here  the  transmitter  was  a  UH-IH  helieopter  equipped  with  4  transmit 
antennas  as  illustrated  in  Figure  8.  The  helieopter  flight  path,  the  multipath  environment,  and 
the  loeations  of  the  reeeive  antennas  are  shown  in  Figure  9.  The  ground  station  eomprised  6- 
ft  diameter  and  4-ft  diameter  traeking  dishes,  labeled  RXl  and  RX2,  respeetively,  loeated  as 
shown  in  Figure  8.  The  ehannel  sounding  signal,  oeeupying  50  MHz  of  bandwidth  and  eentered 
at  1800  MHz  (upper  L-band),  was  used  to  probe  the  ehannel.  The  resulting  impulse  responses 
were  resampled  to  20  Msamples/s  to  produee  a  ehannel  for  10  Mbit/s  SOQPSK-TG  sampled  at 
2  samples/bit.  Additional  details  are  available  in  [10]. 

During  the  flight  line  run,  the  UH-IH  maintained  an  altitude  between  15  and  20  ft.  AGL.  The 
reeeive  antennas  were  positioned  on  60  ft.  towers.  The  hangers  along  the  flight  line  are  also  60 
ft.  tall  at  the  roof  eenters.  Consequently,  line-of-sight  propagation  is  not  possible  when  a  hanger 
is  between  the  UH-IH  and  the  reeeive  antennas.  Consequently  this  data  set  represents  a  rieh 
multipath  propagation  environment. 

To  illustrate  the  behavior  of  p,  we  assign  to  hfin)  the  resampled  version  of  he, 2, lit)  [trans¬ 
mission  from  the  nose  antenna]  and  to  h2in)  the  resampled  version  of  he, 2, Ait)  [transmission 
from  the  tail  antenna].  The  forward-to-aft  separation  between  the  two  transmit  antennas  was  400 
inehes  (approximately  61  wavelengths  at  1800  MHz).  When  the  UH-IH  is  level,  tail  antenna  is 
43  inehes  (approximately  61/2  wavelengths  at  1800  MHz)  higher  than  the  nose  antenna.  During 
flight,  the  helieopter  pitehes  forward  (as  suggested  in  Figure  8)  and  this  plaees  the  tail  antenna 
somewhat  more  than  43  inehes  higher  than  the  nose  antenna. 

The  loeations  of  the  two  transmit  antennas  and  the  reeeive  antenna  produee  impulse  responses 
with  signifieant  (and  similar)  degrees  of  ISI  and  similar  energies,  although  there  will  be  loeations 
where  the  propagation  path  between  one  of  the  two  transmit  antennas  and  the  reeeiver  will  be 
more  attenuated  than  the  other  and  where  the  ISI  on  one  ehannel  may  be  mueh  less  severe  than 
on  the  other.  Thus,  this  set  of  39,300  impulse  response  pairs  represents  a  rieh  mixture  of  varying 
signal-to-noise  ratios  and  ISI. 
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The  results  using  the  natural  normalization  are  summarized  by  the  histograms  shown  in  Fig.  10. 
Curiously,  11,454  of  the  39,300  values  of  p  for  Ei,/Nq  =  10  dB  [Fig.  10  (a)]  are  elose  to  either  0 
or  1.  This  means  approximately  29%  of  the  channel  pairs  prefer  the  use  of  a  single  channel  over 
the  use  of  both  channels.  The  remaining  27,86  channel  pairs  minimize  the  mean  squared  error 
by  applying  some  portion  of  the  available  transmit  power  to  both  antennas.  As  the  signal-to- 
noise  ratio  increases  to  20  dB  [Fig.  10  (b)],  the  number  of  occurrences  where  the  single  antenna 
solution  is  best  falls  to  5,366  or  14%  of  the  channel  pairs.  This  indicates  that  as  signal-to-noise 
ratio  increases,  ISI  plays  an  increasingly  dominant  role  in  determining  the  power  allocation  that 
minimizes  the  mean  squared  error.  The  more  similar  a  pair  of  channels  are  from  the  ISI  point  of 
view,  the  more  p  =  ^/ 2  minimizes  the  mean  squared  error.  This  suggests  that  when  the  signal-to- 
noise  ratio  penalty  is  removed  from  the  weaker  of  the  two  channels,  an  even  stronger  preference 
for  p  =  1/2  should  be  observed. 

That  this  is  true  is  documented  by  the  histograms  in  Figure  11.  In  Figure  II  (a),  710  of  the 
39,300  channel  pairs  (2%)  prefer  only  one  of  the  channels  when  E^/Nq  =  10  dB  and  in  Figure  1 1 
(b),  743  of  the  39,300  channel  pairs  (2%)  prefer  only  one  of  the  channels  when  E^/Nq  =  20 
dB. 

In  summary,  the  results  for  a  multipath-rich  environment  show  that  when  signal-to-noise  ratio 
is  moderate  and  the  channels  are  not  normalized  to  equal  energy,  applying  100%  of  the  available 
power  to  one  of  the  two  transmit  antennas,  is  frequently  the  solution  that  minimizes  the  mean 
squared  error.  As  signal-to-noise  ratio  increases,  or  the  signal-to-noise  ratio  penalty  of  the  weaker 
channel  is  removed,  applying  100%  of  the  available  power  to  one  of  the  two  antennas  is  less 
often  the  solution  that  minimizes  the  mean  squared  error.  In  this  case,  as  long  as  the  channels 
have  similar  degrees  of  ISI,  applying  power  to  both  channels  becomes  the  best  thing  to  do. 
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Antenna  2  (Top) 

fuselage  station  =  85 
center  line  =  0" 
water  line  =  85" 


Antenna  4  (Tail) 

fuselage  station  =  400" 
center  line  =  0" 
water  line  =  65" 


Antenna  1  (Nose  i 

fuselage  station  =  0" 
center  line  =  0" 
water  line  =  22" 


Antenna  3  (Belly) 

fuselage  station  =  100" 
center  line  =  0" 
water  line  =  10" 


RX2 


Fig.  8.  A  diagram  of  the  helicopter-to-ground  multi-antenna  communications  link.  Flere  hc,j,i{t)  is  the  continuous-time  impulse 
response  of  the  channel  between  transmit  antenna  i  and  receive  antenna  j  for  i  G  {1,  3,  4}  and  j  €  {1,2}. 
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Fig.  9.  The  flight  path  for  the  UH-lFl  helicopter  and  the  locations  of  the  two  receive  antennas  for  the  multipath  channel 
sounding  experiments  at  Cairns  Army  Airfield,  Ft.  Rucker,  Alabama.  The  2  receive  antennas  were  located  on  towers  at  the 
bottom  of  the  figure. 
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B.  Edwards  Air  Force  Base,  California  -  Takeoff  and  Landing 

This  data  set  comprises  channel  impulse  responses  collected  at  Edwards  Air  Force  Base, 
California.  Here  the  transmitter  was  a  Beechcraft  C-12  equipped  with  two  transmit  antennas 
as  illustrated  in  Figure  12.  The  antennas  are  separated  by  79.25  inches  (approximately  12 
wavelengths  at  1824.5  MHz)  in  the  forward/aft  direction  and  69.5  inches  (approximately  11 
wavelengths  at  1824.5  MHz)  in  the  up/down  direction.  The  ground  station,  equipped  with  an 
5-meter  parabolic  reflector  autotracking  antenna,  was  located  at  Building  4795  (34°  58'  14.52" 
N,  117°  55'  52.68"  W)  as  shown  in  Figure  13.  The  channel  sounding  signal,  occupying  50  MHz 
of  bandwidth  and  centered  at  1824.5  MHz  (upper  F-band)  was  used  to  probe  the  channel.  The 
resulting  impulse  responses  were  resampled  to  20  Msamples/s  to  produce  a  discrete-time  channel 
for  10  Mbit/s  SOQPSK-TG  sampled  at  2  samples/bit.  See  [20]  for  more  details. 

The  focus  of  this  section  is  the  take-off  and  landing  test  points  on  runway  22F  as  shown  in 
the  Figure  13.  Building  5790  is  approximately  8.75  km  from  each  end  of  runway  22F  and  120 
m  higher  in  elevation.  This  permits  a  least  a  partially  unobstructed  view  of  the  aircraft  when 
it  is  on  either  end  of  the  runway.  The  multipath  environment  comprises  ground  reflections  and 
reflections  off  structures  in  the  line-of-sight  path  between  the  aircraft  and  the  ground  station 
tracking  antenna.  Given  the  narrow  beamwidth  of  the  tracking  antenna  (approximately  2.3°  null- 
to-null),  multipath  interference  tends  to  disappear  once  the  aircraft  achieves  a  modest  altitude 
(approximately  175  m  above  the  runway). 

To  illustrate  the  behavior  of  p  is  for  take-off  and  landing,  hi{n)  is  assigned  to  the  resampled 
version  of  the  upper-antenna-to-ground  propagation  path  and  h2{n)  to  the  resampled  version 
of  the  lower-antenna-to-ground  propagation  path.  The  geometry  is  such  that  the  propagation 
path  between  the  lower  antenna  and  the  receiver  is  partially  obscured  by  the  aircraft  fuselage. 
Consequently,  hi{n)  is  accompanied  by  a  higher  signal-to-noise  ratio  than  h2{n).  The  lower 
antenna  is  further  penalized  by  the  presence  of  a  more  direct  path  to  the  ground-based  reflectors 
that  produce  multipath.  Thus  h2{n)  has  more  ISI  than  hi{n).  This  double  penalty  of  lower 
signal-to-noise  ratio  and  higher  ISI  should  skew  the  optimum  power  allocation  in  favor  of  hi{n) 
(that  is,  p  1). 

The  results,  summarized  in  Figures  14  and  15  for  the  12,564  pairs  of  channel  impulse  responses 
show  that  this  prediction  is  mostly  true.  As  expected,  applying  100%  of  the  available  power 
to  the  upper  antenna  dominates  the  solutions  that  minimize  the  mean  squared  error  for  both 
normalization  modes:  11,800  or  94%  for  Ei,/Nq  =  10  dB  and  11,894  or  95%  for  Eij/Nq  =  20 
dB  with  the  natural  normalization;  11,850  or  94%  for  Eb/No  =  10  dB  and  11,894  or  95%  for 
Eb/Nf)  =  20  dB  with  the  equal-energy  normalization.  The  few  alternatives  to  this  solution  are 
those  where  100%  of  the  available  power  is  allocated  to  the  lower  antenna.  This  is  probably 


26 


Antenna  1  (Upper) 

fuselage  station  =  302" 


center  line  =  10"  (left) 
water  line  =  76" 


Fig.  12.  A  diagram  of  the  airborne  platform  for  the  channel  sounding  experiments  at  Edwards  AFB,  California,  showing  the 
locations  of  the  transmit  antennas. 


explained  by  the  relatively  large  spatial  separation  between  the  two  transmit  antennas:  when 
hi{n)  is  experieneing  signifieant  ISI,  h2{n)  ean  be  relatively  ISI-free.  What’s  more,  when  the 
energy  of  h2{n)  is  signifieantly  below  that  of  hi{n),  the  eontribution  of  the  ISI  assoeiated  with 
h2{n)  to  the  overall  mean  squared  error  is  less.  The  figures  show  that  as  signal-to-noise  ratio 
inereases,  the  seenarios  where  p  =  0  minimizes  mean  squared  error  beeome  less  frequent. 

In  summary,  the  take-off  and  landing  seenario  at  Edwards  AFB,  California,  is  eharaeterized  by 
a  geometrie  arrangement  where  one  of  the  two  ehannels  is  inferior  in  both  signal-to-noise  ratio 
and  ISI.  The  solution  that  minimizes  the  mean  squared  error  is  the  one  where  all  the  available 
power  is  applied  to  the  good  antenna,  and  the  pereentage  of  the  eases  where  this  is  the  solution 
inereases  with  signal-to-noise  ratio. 
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Fig.  13.  Map  of  Edwards  AFB,  California,  showing  the  location  of  the  ground  (receiving)  station  at  Building  4795  (upper  left 
comer  of  the  map)  and  the  locations  of  the  aircraft  during  take-off  and  landing  (bottom  right  corner  of  the  map). 
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Fig.  14.  Optimum  power  allocations  in  the  mean-squared  error  sense  for  the  measured  channel  impulse  responses  from  take-off 
and  landing  on  runway  22L  at  Edwards  AFB,  California,  using  the  natural  normalization:  (top)  Eh/No  =  10  dB;  (bottom) 
Eb/No  =  20  dB. 
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Fig.  15.  Optimum  power  allocations  in  the  mean-squared  error  sense  for  the  measured  channel  impulse  responses  from  take-off 
and  landing  on  runway  22L  at  Edwards  AFB,  California,  using  the  equal-energy  normalization:  (top)  Eb/No  =  10  dB;  (bottom) 
Eb/No  =  20  dB. 
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C.  Edwards  Air  Force  Base,  California  -  Black  Mountain  and  Cords  Road 

These  two  data  sets  comprise  channel  impulse  responses  also  collected  at  Edwards  Air  Force, 
California.  The  airborne  transmitter,  ground  station,  and  center  frequency  described  in  Sec¬ 
tion  V-B  were  used  to  obtain  these  data  sets.  What  is  different  is  the  flight  path.  The  two  flight 
paths,  called  Black  Mountain  and  Cords  Road,  together  with  their  geographic  relationship  to  the 
ground  station  (Building  4795)  are  shown  in  Figure  16.  The  Black  Mountain  flight  path  is  an 
east- west  flight  corridor  from  35°  12'  00"  N,  116°  50'  00'  W  to  35°  12'  00"  N,  117°  56'  00' 
W.  The  Cords  Road  flight  path  is  also  an  east-west  flight  corridor  south  of  the  Black  Mountain 
flight  path  and  extends  from  35°  05'  00"  N,  116°  42'  00'  W  to  35°  05'  00"  N,  118°  22'  00' 
W.  For  both  flight  paths,  the  airborne  transmitter  maintained  an  altitude  of  5,000  feet  AMSF. 
These  flight  paths  are  important  in  that  they  capture  the  behavior  of  low  elevation  angle  “up 
and  away”  flight  paths.  An  important  characteristic  of  this  flight  geometry  is  the  low  elevation 
angle  of  the  receive  antenna.  The  elevation  angle  is  low  due  to  the  large  distances  between  the 
airborne  transmitter  and  the  ground  based  receiver.  (Building  4795  is  approximately  32  km  and 
100  km  from  the  west  and  east  ends  of  the  Black  Mountain  flight  path,  respectively;  29  km  and 
111  km  from  the  west  and  east  ends  of  the  Cords  Road  flight  path,  respectively.)  Because  of  the 
low  elevation  angles,  air-to-ground  propagation  is  prone  to  multipath  interference  from  strong 
“ground  bounces”  [21]. 

As  in  Section  V-B,  hi{n)  represents  the  resampled  version  of  the  upper-antenna-to-ground 
propagation  path  whereas  h2in)  represents  the  resampled  version  of  the  lower-antenna-to-ground 
propagation  path.  Because  of  the  geometry,  the  upper  antenna  is  partially  masked  by  the  aircraft 
fuselage.  Consequently,  the  signal  received  from  the  lower  antenna  is  stronger.  Thus  h2{n)  is 
characterized  by  a  higher  signal-to-noise  ratio  than  hi{n).  But  the  lower  antenna  also  illuminates 
the  ground  (the  source  of  most  of  the  multipath  reflections)  with  a  stronger  signal  than  the  upper 
antenna.  Thus  although  h2in)  has  a  higher  signal-to-noise  ratio  than  hi{n),  h2{n)  also  has  more 
ISI  than  hi{n).  In  this  way,  the  Black  Mountain  and  Cords  Road  flight  paths  present  a  scenario 
involving  a  signal-to-noise  ratio  and  ISI  tradeoff.  How  this  tradeoff  interacts  with  p  to  minimize 
the  achievable  mean-squared  error  is  the  open  question  here. 

The  results  for  the  190,716  pairs  of  channel  impulse  responses  for  the  Black  Mountain  flight 
path  are  summarized  by  the  histograms  in  Figures  17  and  18.  The  results  corresponding  to  the 
natural  normalization  are  shown  in  Figure  17.  Here,  the  ISI  is  such  that  the  equalizer  is  able 
to  reduce  the  residual  ISI  to  the  point  where  signal-to-noise  ratio  is  the  dominant  factor.  The 
end  result  is  that  the  ISI,  signal-to-noise  ratio,  and  power  allocation  tradeoff  becomes  mostly 
a  tradeoff  between  signal-to-noise  ratio  and  power  allocation.  This  tradeoff  is  such  that  the 
mean  squared  error  is  minimized  by  applying  100%  of  the  available  power  to  the  lower  antenna 
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[corresponding  to  h2{n)]  in  the  vast  majority  of  the  cases  (152,246  or  80%  for  Eb/No  =  10  dB, 
130,989  or  69%  for  E^/Nq  =  20  dB).  With  the  equal-energy  normalization,  the  signal-to-noise 
ratio  penalty  associated  with  hi{n)  is  removed  and  the  tradeoff  reduces  to  an  ISI,  signal-to-noise 
ratio  tradeoff.  In  this  case,  summarized  in  Figure  18,  the  situation  reverses.  The  mean  squared 
error  is  minimized  by  applying  100%  of  the  available  power  to  the  upper  antenna  in  the  vast 
majority  of  the  cases  (147,738  or  77%  for  E^/Nq  =  10  dB,  155,786  or  82%  for  E^/Nq  =  20 
dB). 

The  results  for  the  154,305  pairs  of  channel  impulse  responses  for  the  Cords  Road  flight  path 
are  summarized  by  the  histograms  in  Figures  19  and  20.  As  expected,  the  results  are  similar  to 
those  from  the  Black  Mountain  flight  path.  For  the  natural  normalization  (Figure  19),  the  mean 
squared  error  is  minimized  in  the  majority  of  the  cases,  by  applying  all  of  the  available  power 
to  the  lower  antenna  (99,884  or  65%  for  E^/Nq  =  10  dB,  92,566  or  60%  for  E^/Nq  =  20  dB). 
What  is  different  here  is  that  in  those  cases  where  this  is  not  true,  the  preference  is  to  allocate 
100%  to  the  upper  antenna.  This  alternate  preference  becomes  less  pronounced  as  signal-to-noise 
ratio  increases.  For  the  equal-energy  normalization  (Figure  20),  the  preference  reverses  and  the 
mean  squared  error  is  minimized  by  applying  100%  of  the  available  power  to  the  upper  antenna 
in  the  vast  majority  of  the  cases  (116,656  or  76%  for  E^/Nq  =  10  dB,  124,985  or  81%  for 
Eb/No  =  20  dB). 

In  summary,  the  Black  Mountain  and  Cords  Road  flight  paths  are  characterized  by  similar 
geometries.  The  channel  between  the  upper  antenna  and  the  receiver,  hi{n),  is  characterized  by 
lower  signal-to-noise  ratio  and  lower  ISI  than  h2{n),  the  channel  between  the  lower  antenna  and 
the  receiver.  Using  the  natural  normalization,  the  equalizer  is  able  to  remove  ISI  to  the  point 
where  signal-to-noise  ratio  is  the  dominant  factor.  In  this  case  applying  100%  of  the  available 
energy  to  the  lower  antenna  is  the  optimum  approach  in  most  of  the  cases.  Using  the  equal- 
energy  normalization,  differences  in  signal-to-noise  ratio  are  removed  so  that  ISI  becomes  the 
dominant  factor.  Because  the  ISI  on  the  channel  from  the  upper  antenna  is  so  much  less  than 
the  ISI  from  the  lower  antenna,  applying  100%  of  the  available  energy  to  the  upper  antenna  is 
the  optimum. 
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Fig.  16.  Map  of  Edwards  AFB,  California,  showing  the  location  of  the  ground  (receiving)  station  at  Building  4795  and  the  flight  paths  for  the  Black  Mountain  and 
Cords  Road  flight  paths. 
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Fig.  18.  Optimum  power  allocations  in  the  mean-squared  error  sense  for  the  measured  channel  impulse  responses  from  the 
Black  Mountain  flight  path  at  Edwards  AFB,  California,  using  the  equal-energy  normalization:  (top)  Eh/No  =  10  dB;  (bottom) 
Eb/No  =  20  dB. 
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Fig.  20.  Optimum  power  allocations  in  the  mean-squared  error  sense  for  the  measured  channel  impulse  responses  from  the 
Cords  Road  flight  path  at  Edwards  AFB,  California,  using  the  equal-energy  normalization:  (top)  Eb/No  =  10  dB;  (bottom) 
Eb/No  =  20  dB. 
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D.  A  Statistical  Model  for  Triply-Selective  Channels 

It  is  interesting  to  compare  these  results  with  what  might  be  inferred  from  using  a  simple 
statistical  channel  model.  To  do  so,  we  use  a  simple  Gaussian  model  for  each  channel  such  as  that 
used  in  [22] .  In  this  experiment,  channel  1  comprises  1 1  IID  zero-mean  complex- valued  Gaussian 
random  variables  and  channel  2  comprises  18  IID  zero-mean  complex- valued  Gaussian  random 
variables.  These  numbers,  11  and  18,  are  the  average  lengths  of  hi{n)  and  h2{n),  respectively, 
in  the  Cairns  Army  Airfield  data  set  examined  in  Section  V-A.  The  channels  were  normalized 
using  the  equal-energy  normalization  described  above  and  39,300  independent  realizations  were 
produced.  The  results  are  summarized  by  the  histograms  in  Fig.  21.  The  temptation  is  to  think 
of  the  optimum  p  as  a  normally  distributed  random  variable,  but  it  should  be  kept  in  mind  that 
this  is  not  the  case  because  0  <  p  <  1.  For  modest  values  of  E^/Nq  we  observe  that  the  mean 
value  of  the  optimum  p  is  about  0.5.  Given  the  fact  that  p  =  0.5  corresponds  to  traditional 
TR-STBC,  we  see  that  the  simple  statistical  model  suggests  that  traditional  TR-STBC  is  the  best 
on  average.  This  is  in  contrast  to  the  conclusion  drawn  from  the  measured  channel  data,  where 
a  strong  preference  for  transmit  selection  diversity  is  observed. 

As  Eb/ No  increases,  the  optimum  the  optimum  value  of  p  appear  to  increase.  This  is  explained 
as  follows:  as  Eb/No  increases,  the  contribution  to  squared  error  from  additive  noise  is  reduced 
and  ISI  starts  to  become  the  dominant  contributor.  The  optimum  p  >  0.5  means  the  system 
prefers  to  allocate  more  energy  to  channel  1  than  channel  2.  This  makes  sense  because  channel 
1  is  shorter,  and  this  tends  to  contribute  less  residual  ISI  at  the  equalizer  output. 
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VI.  Conclusions 

A  criterion  that  allows  one  to  prediet  when  it  is  better  to  use  transmit  seleetion  diversity  (i.e., 
one  transmit  antenna)  or  the  diversity  aehievable  through  TR-STBC  (i.e.,  two  transmit  antennas) 
was  developed.  The  eriterion  is  the  residual  mean-squared  error  at  the  output  of  an  MMSE 
equalizer.  The  residual  mean-squared  error  was  not  only  a  mathematieally  traetable  quantity,  but 
also  an  exeellent  predietor  of  the  eurious  behavior  illustrated  in  Fig.  1. 

These  eoneepts  were  applied  to  a  set  of  measured  ehannel  impulse  responses  eolleeted  by 
sounding  experiments  in  multipath-prone  loeations  at  Cairns  Army  Airfield,  Ft.  Rueker,  Alabama, 
and  Edwards  Air  Foree  Base,  California.  For  eaeh  pair  of  ehannel  impulse  responses,  the  value  of 
p  that  minimized  the  residual  mean-squared  error  (42)  was  eomputed.  The  eomputed  values  for 
p  were  used  to  form  histograms  to  summarize  the  results.  These  results  illustrate  the  following 
points: 

1)  In  a  2-transmit,  l-reeeive  antenna  system  operating  in  a  frequeney  non-seleetive  fading 
environment,  if  the  two  ehannels  have  unequal  gains,  the  optimum  thing  to  do  is  apply 
all  of  the  available  power  to  the  stronger  ehannel.  That  is,  transmit  seleetion  diversity  is 
optimum.  In  eontrast,  on  a  frequeney  seleetive  fading  ehannel,  the  optimum  approaeh  is  to 
apply  power  to  produee  the  best  trade-off  between  SNR  and  ISI.  That  is,  transmit  seleetion 
diversity  may  not  be  optimum.  The  optimum  value  of  p  assoeiated  with  the  GTR-STBC 
system  deseribed  in  this  paper  identifies  the  best  trade-off  between  SNR  and  ISI.  There 
are  some  ehannel  pairs  for  whieh  p  =  0  or  1  is  the  optimum  (transmit  diversity  ease)  and 
some  ehannel  pairs  for  whieh  p  =  ^/2  (traditional  TR-STBC).  But  there  are  many  ehannel 
pairs  for  whieh  neither  of  these  is  optimum. 

2)  On  our  measured  ehannels,  transmit  seleetion  diversity  was  more  eommon  than  traditional 
TR-STBC. 

3)  On  a  statistieal  ehannel,  sueh  as  the  one  used  in  [22],  traditional  TR-STBC  is  the  best 
thing  to  do  on  average.  This  is  in  eontrast  to  the  results  from  the  measured  ehannels. 
Consequently,  the  optimum  power  alloeation  in  a  real  setting  is  not  predieted  well  by 
simple  statistieal  ehannel  models. 

The  expression  for  the  residual  mean-squared  error  (42)  was  developed  by  formulating  a 
generalization  of  TR-STBC  that  ineluded  not  only  transmit  seleetion  diversity  and  traditional 
TR-STBC,  but  also  a  generalization  that  permitted  the  transmitter  to  alloeate  unequal  power 
to  eaeh  transmit  antenna.  This  generalization,  ealled  generalized  TR-STBC  (GTR-STBC),  ean 
be  used  as  the  basis  of  transmit  diversity  system  where  the  transmitter  has  partial  ehannel 
state  information  (in  the  form  of  p).  This  transmit  diversity  system  is  apropos  to  single  earrier 
modulations  operating  in  situations  with  peak-to-average  power  ratio  eonstraints. 
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